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We have recently proposed a dynamical mechanism that may realize a flat four-dimensional space 
time as a brane in type IIB superstring theory. A crucial role is played by the phase of the chiral 
fermion integral associated with the IKKT Matrix Theory, which is conjectured to be a nonper- 
turbative definition of type IIB superstring theory. We demonstrate our mechanism by studying a 
simplified model, in which we find that a lower-dimensional brane indeed appears dynamically. We 
also comment on some implications of our mechanism on model building of the brane world. 



PACS numbers ll.25.-w; 11.25.Sq 

Introduction. — The idea that our four-dimensional 
space time is realized as a brane in a non- compact higher- 
dimensional space time has recently attracted much at- 
tention. Through many works during the last few years, 
it is expected to provide natural resolutions to many long- 
standing problems in the Standard Model. The hierarchy 
problem is transmuted into a geometrical one [JjJ, and 
it was further argued that the exponential dependence 
of the "warp" factor in the extra directions reduces the 
problem to a fine-tuning of order 50 J3 . The cosmological 
constant problem may also be resolved in such a setup || . 
It has been argued that any nonzero higher-dimensional 
cosmological constant is absorbed into the warp factor, 
and that the four-dimensional cosmological constant is 
automatically tuned to zero (or to a very small number). 
A possible obstruction to the idea (as opposed to a more 
conventional idea using Kaluza-Klein compactifications) 
is that gravity may propagate in higher dimensions and 
thereby contradicts the 4D Newton's law observed in the 
low energy scale. However, the particular AdS-type back- 
ground metric that arises naturally in such a setup allows 
a normalizable zero mode of the graviton bound to the 
brane Q . A small correction to the 4D Newton's law due 
to the continuum spectrum of massive modes is argued to 
be small enough to be compatible with the experimental 
bound. All these attractive features of the idea lead us 
to hope that there is a natural string-theory realization 
of the brane world scenario. 

In Ref. [|, we have proposed a dynamical mechanism 
which may realize a flat four-dimensional space time as 
a brane in type IIB superstring theory. Obviously, such 
a mechanism should inevitably be of nonperturbative na- 
ture. Indeed, our mechanism was based on the IKKT 
version Q of the Matrix Theory [Q , namely the IIB ma- 
trix model, which is conjectured to be a nonperturbative 
definition of type IIB superstring theory. The model is a 
supersymmetric matrix model composed by ten N x N 
hermitian bosonic matrices and sixteen N x N hermi- 
tian fermionic matrices. The space time is represented 
by the eigenvalues of the bosonic matrices. The model 
has manifest ten-dimensional Lorentz invariance, where 
the bosonic and fermionic matrix elements transform as a 



vector and a Majorana-Weyl spinor, respectively. The in- 
tegral over the fermionic matrices yields a pfaffian which 
is complex in general. This poses a technical difficulty 
known as the 'complex action' problem in studying the 
IIB matrix model by Monte Carlo simulation. Monte 
Carlo studies incorporating only the modulus of the pfaf- 
fian (and omitting the phase by hand) showed that the 
space-time becomes isotropic in ten dimensions in the 
large- N limit || |J. This result suggests that the phase 
of the pfaffian must play a crucial role if a brane world 
naturally arises in the type IIB superstring theory. The 
effect of the phase is to favour configurations for which 
the phase becomes stationary. Such an effect has been 
studied within a saddle-point approximation and found 
to enhance lower-dimensional brane-like configurations 
considerably 

In this Letter, we demonstrate our mechanism more 
explicitly by studying a simplified model using Monte 
Carlo simulation. In this case, we find that the dominant 
saddle-points are given by configurations with only three- 
dimensional extent. 

The mechanism. — The IIB matrix model fl is for- 
mally a zero- volume limit of D = 10, J\f = 1, pure super 
Yang-Mills theory. The partition function of the IIB ma- 
trix model (and its generalizations to D = 4 and D = 6) 
can be written as 

Z im = JdAe- s * Z f [A] , (1) 

where S h = -Tr (L4 M , A u } 2 )/4, and T[A] = -\nZ f [A} 
represents the effective action induced by integration over 
the fermionic matrices. The dynamical variables A^ (fx — 
1, • • • , D) are D bosonic N x N traceless hermitian ma- 
trices. Expanding A^ as A^ = X^i" 1 A®t a in terms of 
the generators t a (a = 1, . . . , N 2 - 1) of SU(iV), the inte- 
gration measure dA is given as dA = Yla=i 11,^= i" 1 dA^. 
The generators t a are normalized as Tr (t a t h ) = 25 a b- 

The fermion integral Zf[A] is complex in general for 
D = 10, N > 4 and for D = 6, N > 3 @ jn|. We 
restrict ourselves to these cases in what follows. In the 
D = 10 case, the fermion integral Zf[A] is given by the 
pfaffian FfM, where M is a 16(A^ 2 - 1) x 16(N 2 - 1) 
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complex antisymmetric matrix defined by 
M aa ,bf3 = Tr (t a (CT^) af3 [A tt ,t b ]) 



(2) 



regarding each of (aa) and (6/3) as a single index. Here, 
{[i — 1, . . . , 10) are ten-dimensional Weyl-projected 
16 x 16 gamma matrices, satisfying CT^C^ — (r M ) T with 
C = C T being the unitary charge conjugation matrix. 
Similarly in the D = 6 case, the fermion integral Zf[A] 
is given by the determinant det where is a 

4(7V 2 — 1) x 4(A 2 — 1) complex matrix defined by 



M { Z b ^Tr{t a (r^U[A,,t b } 



(3) 



regarding each of (aa) and (6/3) as a single index. Here, 
I/j (a 1 = 1, - - ■ , 6) are six-dimensional Weyl-projected 
4x4 gamma matrices. 

Since the fermion integral Zf[A] is a complex quan- 
tity for the cases under consideration, let us write it as 
Z f [A] = exp(-rM - ffW). In Ref. §, the effect of 
the phase in the path integral has been studied 
using a saddle-point approximation, whose validity has 
been also discussed. The saddle-point equation for rM 
is given by 
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It is useful to introduce the following classification of 
"brane" configurations 
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where njj (i — 1, ■ ■ ■ , D — d) are (D — d) linearly inde- 
pendent D-dimensional real vectors. Namely, fid repre- 
sents a set of configurations with less than d-dimensional 
extent. Note that f2i C £1% C • • • C flu, where flo is no- 
thing but the whole configuration space of the model. In 
Ref. p] we proved that all configurations in £Id-2 are so- 
lutions to the saddle-point equation (^). Assuming that 
the configurations in are the dominant saddle-point 

configurations, we still have to integrate over those con- 
figurations to determine the actual dimensionality of the 
space-time. In fact, the gaussian fluctuation of the phase 
around the saddle-points gives a huge enhancement 
to the brane configurations with lower dimensionality, 
and this enhancement cancels exactly the entropical bar- 
rier against having such configurations. In the D = 10 
case, this provides a dynamical mechanism for the possi- 
ble appearance of four-dimensional space time as a brane 
in ten-dimensional space time. 

A simplified model. — In order to investigate how our 
mechanism works, we consider a simplified model descri- 
bing an integration over the saddle-point configurations. 
Specifically, we consider the integral 



Z = dA 



3 -/3F[A]- 7 /[A] 



where the functions F[A] and f[A] are defined as 
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Since the function F[A] vanishes if and only if the con- 
figuration {An} satisfies the saddle-point equation (Q), 
the integral m) is dominated by the saddle-point config- 
urations in the large- (3 limit. The function f[A] makes 
the integral (||) convergent as long as the coefficient 7 
is fixed to be a real positive number. In fact, the para- 
meter 7 can be absorbed by an appropriate rescaling of 
A^ and (3. Therefore we take 7 = 1/2 in what follows 
without loss of generality. Note also that the simpli- 
fied model (^|) is invariant under a SO(D) transformation 
A,j, h-> A.n V A v , where A M „ G SO(D), and a SU(A) trans- 
formation A^ 1— > gA^g' , where g € SU(iV), which are the 
symmetries of the original model ([!]). 

Using the invariance of the partition function 
the change of variables A® 
exact relation 
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XAf., one can obtain an 



(9) 



where the ensemble average ( ■ )p is defined with the par- 
tition function (||). Assuming that {f)p goes to a con- 
stant c for (3 — > 00, we obtain the asymptotic behavior of 
(F)f3 for large /? as 
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where the coefficient b is given as b = c/2 — D{N 2 — 
l)/2. This confirms the above claim that the integral (||) 
is dominated by the saddle-point configurations in the 
large-/? limit. 

A quantity which fully characterizes the dimensionality 
d of a given configuration can be given by the moment 
of inertia tensor T defined by the D x D real symmetric 
matrix M] 



T — 
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A configuration {A^} belongs to fi^, if and only if the 
number of zero eigenvalues of the matrix T is more than 
or equal to (D — d). Let us denote the eigenvalues of the 
matrix T as (i — 1, • • • , D), where Ai > A2 > • • • > 
Ad > 0. We can determine the dimensionality of the 
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dominant saddle-point configurations from the ensemble 
average of the eigenvalues (Aj)/3 in the (3 — > oo limit. 

We address this issue by performing Monte Carlo 
simulation using a Metropolis algorithm. We create a 
trial configuration {A'^} by replacing an element of 
the previous configuration {A^} with a new one gener- 
ated with the probability distribution exp[— i(A^) 2 ]. 
The trial configuration is accepted with the probability 
min(l,exp(-/3AF)), where AF = F[A'] - F[A]. This 
procedure is repeated for all the elements of the con- 
figuration. The computational effort required for the 
above algorithm is of order 0(iV 8 ) per one sweep, which 
is much larger than that for the simulation encountered 
12 1 . Due to this, results with high statistics 



of the simplified model (|6J) has only three-dimensional 
extent. Preliminary results for D = 10, N = 4 suggest 
that this is the case also for D = 10. 
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are obtained only for the case of D = 6 and N = 3 (we 
have made 192,000 accepted updates for each j3 < 384 
and 768,000 for /3 = 512). 




FIG. 1. In the upper part, we plot the left hand side of 
eq. Q against f3 for D — 6 and N = 3. The dashed line rep- 
resent the exact result ^-{N 2 — 1) = 24. In the lower part, the 
function {^[^4])^ is plotted against (3 in a log- log scale. The 
straight line represents a fit to the predicted large-/? behavior 
(F[A])p ~ b/p. 

Results. — In the upper part of Fig. [j] we plot the left 
hand side of (^|) , which demonstrates the validity of our 
simulation. In the lower part of Fig. [l]we plot the average 
{F[A])p against (3 in a log-log scale. The straight line 
represents the fit of the data for (3 > 16 to the predicted 
large-/? behavior @ with b = 1.7(1). 

Fig. H shows the six eigenvalues of the moment of in- 
ertia tensor T as a function of (3. We find that the three 
smallest eigenvalues {i = 4,5,6) are monotonously 

decreasing with a pronounced power law behavior. Fit- 
ting the data for (3 > 16 to the power law behavior, we 
extract the powers -0.040(4), -0.199(2) and -0.450(6), 
respectively. Similarly, the powers for TV = 4 are ex- 
tracted to be -0.11(1), -0.26(2) and -0.36(2). Thus we 
conclude that the dominant saddle-point configurations 




FIG. 2. The six eigenvalues (A^)^ of the moment of inertia 
tensor T are plotted against in a log-log scale for D = 6 and 
N = 3. The straight lines for the smallest three eigenvalues 
represent the fits to the power law behavior. 



Discussion. — The results presented in the previous 
section shows clearly that the stationarity of the phase 
indeed enhances lower-dimensional brane configura- 
tions considerably, thus demonstrating our mechanism. 
In particular, while the existence of saddle-point config- 
urations other than the configurations in £Id-2 is not 
excluded, our results suggest that such configurations, 
even if they exist, can safely be neglected on statisti- 
cal grounds. Given this observation, the reason why 
we obtain the dimensionality 'three' from the dominant 
saddle-point configurations of the simplified model (|^) 
can be understood analytically. As has been done in 
Ref. JB| for the IIB matrix model, we can rewrite the 
[3 — > co limit of the simplified model (^) as an integral 
over the configurations in r2 J o_2 ■ The gaussian fluctua- 
tion of the phase should be taken into account by 
the corresponding Hesse matrix, which is, in the present 
case, just the square of the one for the IIB matrix model 
(or its D = 6 version). Due to this, the gaussian fluc- 
tuation enhances lower-dimensional brane configurations 
even more strongly than in the IIB matrix model and 
overwhelms the entropical barrier against having those 
configurations. 

The fact that the lowering of the dimensionality stops 
at three instead of continuing further down can be under- 
stood as follows. We first note that the fermion integral 
Zf[A] vanishes for configurations in ^2 ||. Therefore, 
the phase is actually ill-defined for configurations in 
r^2- Still, we can consider configurations with A\, A2 
being generic and Ai (i — 3, • ■ • , D) being of order e. 
One can easily see that the function F[A] is diverging as 
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e~ 2 for e — > 0. Therefore, the /3 term in (g) suppresses 
configurations in Q2 strongly. (In other words, configura- 
tions in are not saddle-point configurations, although 

^2 C Sl D -2-) 

As is clearly shown in the present work, the enhance- 
ment occurs precisely when the space-time becomes a 
flat lower-dimensional hyperplane, which we consider as 
a very attractive feature of our mechanism. Namely, our 
mechanism has a built-in structure in which the brane 
that appears as a result of the nonpcrturbative string dy- 
namics is very likely to be flat. According to our mecha- 
nism, scenarios with two separated branes (i.e. our world 
and the so-called 'Planck' brane as in Ref. ||), their ex- 
tensions to many brane s p3[ , and scenarios with mutu- 
ally intersecting branes H] seem to be unnatural. 

Let us also comment on a connection of our mechanism 
to the brane world scenario. In Ref. (ll|, the IIB matrix 
model is expanded around a D3-brane configuration per- 
turbatively and four-dimensional noncommutative Yang- 
Mills theory has been obtained Jl(| . The (perturbatively 
stable) theory, which is obtained in this way from the 
IIB matrix model, has been recently identified jfl| with 
a type IIB superstring theory in AdSs x S 5 with an in- 
finite -B-field background. Remarkably the metric that 
appears in the corresponding supergravity solution takes 
the form of Randall-Sundrum's type and thereby al- 
lows for a four-dimensional Newton's law. We expect 
that brane configurations similar to the D3-brane con- 
figuration considered above as a background in the IIB 
matrix model should appear dynamically as a result of 
our mechanism. Thus our mechanism is rather directly 
related to the brane world scenario. 

In the IIB matrix model, the enhancement and the 
entropical barrier are exactly balanced and the actual di- 
mensionality of the brane is expected to be determined 
as a result of large N dynamics. In this regard, we recall 
that the low-energy effective theory of the IIB matrix 
model has been shown to be described by a branched- 
polymer like system in Ref. [fl9|| . There it was further 
argued that a typical double-tree structure that appears 
in the effective theory might cause a collapse of the con- 
figuration to a lower-dimensional manifold. Whether the 
actual dimensionality of the brane turns out to be four 
or not can be investigated directly by performing the in- 
tegration over the saddle-point configurations as formu- 
lated in Ref. pj. We hope that Monte Carlo techniques 
developed in Ref. J^[ will enable us to address such an 
issue in near future. 
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